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Abstract 

We examine the accuracy of estimation of parameters of the gravitational-wave signals from 
spinning neutron stars that can be achieved from observations by Earth-based laser interferometers. 
We consider a model of the signal consisting of two narrowband components and including both phase 
and amplitude modulation. We calculate approximate values of the rms errors of the parameter 
estimators using the Fisher information matrix. We carry out extensive Monte Carlo simulations 
and obtain cumulative distribution functions of rms errors of astrophysically interesting parameters: 
amplitude of the signal, wobble angle, position of the source in the sky, frequency, and spindown 
qq ' coefficients. We consider both all-sky searches and directed searches. We also examine the possibility 

, of determination of neutron star proper motion. We perform simulations for all laser-interferometric 

detectors that are currently under construction and for several possible lengths of the observation time 
and sizes of the parameter space. We find that observations of continuous gravitational-wave signals 
from neutron stars by laser-interferometric detectors will provide a very accurate information about 
their astrophysical properties. We derive several simplified models of the signal that can be used in 
the theoretical investigations of the data analysis schemes independently of the physical mechanisms 
generating the gravitational-wave signal. 



1 Introduction 



Detection of gravitational waves from spinning neutron stars by currently constructed long-arm laser 
interferometers |], ^ 0, EJ is expected to provide a wealth of astrophysical information about these 
objects. In the first paper of this series || (hereafter Paper I) we introduced a model of the gravitational- 
wave signal from a spinning neutron star. The response of a laser interferometer to such a signal will 
carry information about the neutron star period of rotation and its evolution, amplitude of the signal 
and its polarization. Moreover even the response of a single detector will contain the full information 
about the position of the source in the sky. Also higher order terms in the phase contain velocity of 
the source and its distance. In this second paper of our series we investigate in detail how accurately 
these astrophysically interesting parameters can be determined from observations by laser-interferometric 
detectors. We assume that optimal data analysis methods presented in Paper I are used. We use a 
standard tool, the Fisher information matrix, to assess the rms errors of our estimation method. 

We find impressive potential accuracies achievable with gravitational-wave detectors that can make 
them into astronomical laboratories competitive with optical and radio observatories. For all-sky searches 
and 120 days of observation time the initial laser-interferometric detectors will be able to locate the 
strongest fast and young neutron star gravitational-wave emitters at a distance of 1 kpc with accuracy 
of order 10 -7 sr and estimate the frequency and the four spindown parameters of the wave with relative 
accuracies of order 10~ 2 ( 4 ~ fe ) j where k = 0, ... ,4 denotes the fcth spindown (k = corresponds to the 
frequency of the wave). Amplitude can be estimated within a few percent and wobble angle within a 
few percent of a radian. For directed searches the relative rms errors of estimation of frequency and 
spindowns will decrease with respect to the all-sky searches by a factor of 10 for the fourth spindown, 
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10 2 for the third and the second spindowns, 10 3 for the first spindown, and 10 4 for the frequency. With 
advanced detectors a reasonable accuracy of estimation of proper motions of neutron stars with all-sky 
searches will be achievable for nearby neutron stars (at distances of 40 pc or so). 

The plan of the paper is as follows. In Sec. 2 we review the model of the gravitational-wave signal 
introduced in Paper I. We introduce eight special cases of the model depending on whether observation 
time is long or short and neutron star is young or old and it is spinning fast or slow. In Sec. 3 we 
present calculations of the rms errors that can be achieved for our models assuming that we use maximum 
likelihood detection as data analysis method. To estimate the rms errors we use Fisher information matrix. 
In Sec. 4 we explore the possibility of estimating proper motion of neutron stars from gravitational- wave 
observations. In Sec. 5 we introduce a number of simplified signal models that enable us to assess the rms 
errors of the phase parameters independently of the physical mechanisms generating the gravitational 
radiation from a spinning neutron star. In Sec. 6 using a simplified model we investigate the dependence 
of the rms errors of the parameters on the observational parameters: duration of observation, initial 
moment of observation and the latitutde of the detector. 

Several additional topics are investigated in the appendices. In Appendix A we interpret our 1/4 of 
a cycle criterion in terms of the signal-to-noise ratio loss. In Appendix B we give details of number of 
cycles calculations for our special models. In Appendix C we study the polynomial phase model and the 
effect of the choice of the initial time in the parameter estimation problem. 



2 A model of the two-component gravitational-wave signal 

In Sec. II of Paper I we have introduced the following two-component model of the gravitational-wave 
signal h from a spinning neutron star: 



h(t) = fti(t) + h 2 (t), (1) 



with 



hx(t) = F + (t)h 1+ (t) + F x (t)h lx {t), h 2 {t) = F+{t)h 2+ {t) + F x (t)h 2x (t), (2) 
h 1+ {t) = ■|ft, sin26»sin2/,cos*(i), h 2+ (t) = \h Q sin 2 0(1 + cos 2 l) cos 2*(f), (3) 
hi x (t) = \h sin 29 sin t sin ^ (t) , h 2 x (t) = h Q sin 2 9 cos u sin 2^ (t) . (4) 

The model of the signal given above represents the quadrupole gravitational wave that is emitted by a 
freely precessing axisymmetric star. The angle 9, called the wobble angle, is the angle between the total 
angular momentum vector of the star and the star's axis of symmetry and l is the angle between the 
total angular momentum vector of the star and the direction from the star to the Earth. For the case 
when 9 = ir/2 (then the component hi vanishes) the model also gives the quadrupole wave from a triaxial 
ellipsoid rotating about a principal axis. In both cases the amplitude h a is given by 

k o = —A — > ( 5 ) 

where / is the moment of inertia with respect to the rotation axis and r a is the distance to the star. For 
precessing axisymmetric star f Q is the sum of the frequency of rotation of the star and the frequency of 
precession and e is the poloidal ellipticity of the star whereas for triaxial ellipsoid f is the frequency of 
rotation of the star and e is ellipticity of the star defined by e := (I\ — I 2 )/I, where 1\ and I 2 are the 
moments of inertia of the star with respect to the principal axes orthogonal to the rotation axis. 
Inserting numerical values into Eq. (|^) gives 

^- 4 - 23xl0 ^(iik) 2 ' ( 6 ) 

where 

, = _J I 1 kpc 

° ' 10- 5 10 45 gem 2 r a ' K ' 

In all numerical simulations presented in the next sections except for the case when the effect of the 
proper motion was studied we had assumed the numerical value of the parameter d Q to be 1. The most 
uncertain parameter in Eq. ([?]) is the ellipticity e, the value of e ~ 10~ 5 is thought to be an upper limit 
and typical values may be significantly less || . 

The beam-pattern functions F + and F x from Eq. (|^) are given by 

F+(t) = sinC[a(t)cos2^ + 6(t)sin2^], (8) 
F x (t) = sin ([b(t) cos 2ip - a(i) sin 2^] , (9) 
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where 

a(t) = — sin 27(3 — cos2A)(3 — cos2<5) cos[2(a — <fr r — Q r t)] 

— — cos27sinA(3 — cos 25) sin[2(a — <j) r — Q r t)] 
+ — sin 27 sin 2A sin 25 cos [a — <f> r — fi r t] 

— — cos 27 cos A sin 25 sin [a — 4> r — fi r t] 
3 

+— sin 27 cos 2 A cos 2 5, (10) 
b(t) — cos27sin Asin<5cos[2(a — <f> r — Q, r t)\ 

+— sin 27(3 — cos2A) sin5sin[2(a — <\> r — f2 r t)] 
+ cos 27 cos A cos 5 cos [a — <j> T — f2 r t] 

+— sin 27 sin 2A cos 5 sin [a — <fi r — f2 r t]. (11) 

The angles a, 5, and "0 are respectively right ascension, declination of the gravitational-wave source, and 
polarization angle (these angles determine the orientation of the wave reference frame with respect to 
the celestial sphere reference frame). The angle A is the latitude of the detector's site, 7 determines the 
orientation of the detector's arms with respect to local geographical directions and C is the angle between 
the interferometer arms. Q r is the rotational angular velocity of the Earth, and <j) r is a deterministic 
phase which defines the position of the Earth in its diurnal motion at t = 0. The derivation of the 
formulae (§|)-([ll]) is given in Sec. II A of Paper I. 
The phase \& of Eqs. ||) and (§|) is given by 

*(*) = $„ + 2^5>t— — + —no • r ES (t) £ f - + -n • r E (*) £ /„-, (12) 

fc=0 1 '' fc=0 ' k=0 

where tes is the vector joining the solar system barycenter (SSB) with the center of the Earth and te joins 

(*) 

the center of the Earth with the detector, f a is the kth time derivative of the instantaneous frequency at 
the SSB evaluated at t = 0, no is the constant unit vector in the direction of the star in the SSB reference 
frame. We have neglected the relativistic effects. To derive the above model we assumed that in the rest 
frame of the neutron star the gravitational-wave frequency can be expanded in a power series. For the 
detailed derivation of the phase model ([l2] ) see Sec. II B and Appendix A of Paper I. 

Neglecting the eccentricity of the Earth's orbit and the motion of the Earth around the Earth-Moon 
barycenter the scalar products n • r E and n • r E s of Eq. ( [l2| ) can be written as (see Sec. II B of Paper I) 

no • te = Re [sin A sin 5 + cos A cos 5 cos (a — 4> r — Vt r i)] , (13) 
no • tes = Res [cos a cos 5 cos (<j> + Q t) + (cos e sin a cos <5 + sine sin (5) sin (0 O + Q t)] , (14) 

where Res — 1 AU is the mean distance from the Earth's center to the SSB, Re is the mean radius of 
the Earth, Q a is the mean orbital angular velocity of the Earth, and <p is a deterministic phase which 
defines the position of the Earth in its orbital motion at t — 0, e is the angle between ecliptic and the 
Earth's equator. 

(*) 

The extremal values of the spindown parameters f a arc estimated from the formula 
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/„- (-l) fc fc!^, fe = l,2,..., (15) 

where r we call the spindown age of the star. 

We choose the number of terms in the power series on the right-hand side of Eq. ([l2]) according to 
the following criterion: we exclude an effect from the model of the signal in the case when it contributes 
less than 1/4 of a cycle to the phase of the signal during the observation time. An interpretation of the 
1/4 of a cycle criterion in terms of the signal-to- noise ratio loss is given in Appendix A. If we restrict to 
observation times T < 120 days, frequencies f a < 1000 Hz, and spindown ages r > 40 years, the phase 
model ( |l2| ) meets the criterion for an appropriate choice of the numbers s l5 S2, and S3. Also the effect 
of the star proper motion in the phase is negligible if we assume that the star moves w.r.t. the SSB not 
faster than 10 3 km/s and its distance to the Earth r a > 1 kpc. However for nearby neutron stars at 



distances r Q > 40 pc the proper motion contribution to the phase (12) can exceed 1/4 of a cycle. We 
discuss this issue at the end of the present section. 
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Table 1: The number of spindown terms needed in various contributions to the phase of the signal 
depending on the type of population of neutron stars searched for. The number si refers to the spindown 
contribution, S2 refers to the Earth orbital motion, and S3 refers to the Earth diurnal motion contribution 
[cf. Eq. H)]. 



We consider several models depending on the range of the gravitational- wave frequency f g [f g = f a for 
the hi component and f g = 2f for the h% component of the signal (Q)], spindown age r, and observation 
time T a . Following Brady et al. || we say that neutron star is slow if f g ~ 200 Hz and fast if f g ~ 1 kHz. 
We define a neutron star to be old if r ~ 10 3 yr and young if r ~ 40 yr. We consider observation time to 
be short if T ~ 7 days and long if T Q ~ 120 days. Accordingly we have eight possible models. In Table 1 
we present the number of terms needed in the power series of various contributions to the phase in order 
to meet the above criterion. The details of the calculation are given in Appendix B. 

In our calculations of the rms errors of the estimators of the parameters we have also studied the case 
of directed searches. By directed searches we mean the search for the signal from a spinning neutron star 
which position in the sky is known. This reduces the dimension of the parameter space by 2 since the 
right ascension and the declination of the source are known. We do not assume that the gravitational 
wave frequency is known. Even if we are searching for a gravitational wave from a known pulsar its 
frequency may not coincide with the frequency of the electromagnetic radiation or be exactly twice that 
frequency For directed searches the values of the parameters s\, S2, and S3 given in Table 1 remain in 
general the same however for specific positions in the sky they may be less. 

The numbers of spindowns that must be included in the phase model given in Table 1 agree with the 
corresponding numbers of Brady et al. [0] obtained by a different criterion for the cases they considered 
(see Figure 6 of their paper). 

We have also considered the effect of the proper motion of the neutron star on the phase of the signal 
assuming that it moves uniformly with respect to the SSB reference frame. We have found that for the 
observation time T — 120 days and the extreme case of a neutron star at a distance r = 40 pc, transverse 
velocity |v ns j_| = 10 3 km/s (where v n s_L is the component of the star's velocity v ns perpendicular to the 
vector no), gravitational-wave frequency f g = 2f a = 1 kHz, and spindown age r = 40 years, proper 
motion contributes only ~4 cycles to the phase of the signal. This dominant proper motion contribution 
f pm to the phase ( |T^ ) is given by (cf. Eq. (A20) in Appendix A of Paper I) 

*pm(*) = Jot. (16) 

c r a 

The ratio v ns j_/r determines the proper motion of the star. We shall study the accuracy of estimation 
of the proper motion in Sec. 4. 



3 Estimation of parameters of the two-component signal 

In Sec. Ill of Paper I we have presented an optimum data analysis method: maximum likelihood detection. 
This method consists of maximizing the likelihood function A with respect to the parameters of the signal. 
If the maximum of A exceeds a certain threshold calculated from the false alarm probability that we can 
afford we say that the signal is detected. The values of the parameters that maximize A are said to be 
the maximum likelihood estimators of the parameters of the signal. Let us collect the signal parameters 
into the vector 6 = (0\, . . . ,6 n ). The covariance matrix C of the estimators of the parameters 6 is 
approximately given by the inverse of the Fisher information matrix Y: C = T . This approximation is 
the better the higher the signal-to-noise ratio. For a signal h buried in a Gaussian noise the components 
of the Fisher matrix Y are given by 

/ £lh PlU \ 

(17) 
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where the scalar product ( ■ | • ) is defined by 

(x\y):=mj™^^±df, (18) 

where ~ denotes the Fourier transform, * is complex conjugation, and Sh is the one-sided spectral density 
of the detector's noise. 

For a signal h = hi + h 2 that consists of two narrowband components around the frequencies f {hi 
component) and 2f a {h 2 component), assuming that over the bandwidth of the signal the noise spectral 
density Sh is nearly constant and equal to Sh{fo) or Sh(2f ), the Fisher matrix L is approximately given 

by 

r To/2 dhi dh x , 2 f To/2 dh 2 dh 2 



Tij ~ s h (f ) J_ To/2 de, de j dt + s h (2f a ) J_ To/2 de, oe dt > (19) 

where T is the observation time and the observation interval is [— T /2, T /2). 

For the two-component signal defined by Eqs. (0)-(§) we can extract explicitely the dependence of 
the elements of the covariance matrix C on the amplitude parameter h Q . Let us write the signal ([l]) in 
the form 

h(t ; e) = h [ X i(t;C) + X2(t;C)}, (20) 

where 6 = (h Q , £), i.e. £ is the vector of all the signal parameters with the exception of h . Let us note 

(i) W p—. 

that £ = (a, S, tp, 1, 8, $oi J01 fo>- ■ ■ > fo), so the signal (ECj) depends on 8 + s parameters, where s is the 

number of the spindowns one wants to include. Substituting Eq. ( p0| ) into (|l^) one can show that the 
components of the covariance matrix C = T^ 1 can be written using the matrix notation as (superscript 
T denotes matrix transposition): 

C hoho = K, (21) 
C KC £* -±KA'(A")-\ (22) 

fl 

C cc - ^[(A'O^ + ^A'O-^A'fA'CA")- 1 ], (23) 



where 



K := [A - A'(A") _1 (A') T ] \ (24) 

9 fTo/2 9 fTo/ 2 

A ;= TTn / ixi(t;C)} 2 dt + —— [x 2 (t;C)] 2 ^ (25) 

dh\Jo) J-To/2 dh{4Jo) J-To/2 

A" 2 f T °/ 2 d X i(t;C)d X i(t;C) , , 2 Z^ 2 % 2 (t; Q d X2 (t; Q 

u ' S h (f )J_ To/2 dQ dQ S h (2f )J_ To/2 dQ dQ ' 1 1 



The optimal signal-to-noise ratio d := \J (h\h) for the signal ( |20|) can be computed from the formula 
(cf. Sec. Ill C of Paper I): 

r,u2 rTo/2 9 , 2 f To/2 

[Xi(t-Xfdt + —^- [X2{t-X)f dt. (28) 



Sh(fo) J -To/2 ' Sh(2fo) J -To/2 

We have used the formulae (^l])-(|27]) to calculate numerically covariance matrices for different phase 
models. Because the components of the covariance matrices depend in a complicated way on the trigono- 
metric functions of the angles a, 8, ip, 1, and 9 we have resorted to the Monte Carlo simulations. For 
each simulation of a covariance matrix we have generated 1000 sets of the angles {a, 5, tp, 1, 9} according 
to the probability measure 

da x d sin 5 x dtp x d cos t x d9 

defined on the parameter space 

{a e [0,27r),sin£ € [-1,1], [0,27r),cost e [-1,1], 6 € [0,tt]}. 

In the simulations we have assumed the following models of the noise spectral densities Sh in the 
individual detectors. The noise curves for the VIRGO and the initial/advanced LIGO detectors are 
adopted from 0, and the noise curve for the TAMA300 detector is taken from ||, ||. Wideband and 
narrowband versions of the GEO600 detector noise are based on [91 . 
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We present results of the simulations by plotting the cumulative distribution functions for the square 
roots of the diagonal components of the covariance matrix which approximate the rms errors of the 
estimators of the parameters. For some cases we also give tables with quartiles of the distributions po| . 

We study the rms errors (denoted by a) of the following astrophysically important parameters: amplitude 

(*) 

h a , wobble angle 9, position of the source on the sky, frequency f a , and spindown parameters f a . The 
accuracy of the position on the sky is measured by the solid angle Afi corresponding to the ellipse of 
semiaxes <r(5) and cr(a) and is given by 

Afl = t: cos 5a(5)a(a). (29) 

(k) 

To assess the relative rms errors of the spindown parameters f a we have assumed in the simulations the 
extremal values of these parameters given by formula (fL5|) . 

In Figure 1 we have given cumulative distributions for the rms errors of the parameters for all the 
laser-interferometric detectors under construction. For the GEO600 detector we have also considered the 
narrowband configuration tuned to 1 kHz with bandwidth of 30 Hz. We have taken long observation time 
(T = 120 days) and we have assumed the neutron star to be fast (/ Q = 500 Hz) and young (r = 40 yr). 
In addition in Table 2 we have given the quartiles of the distributions plotted in Figure 1. We observe 
that the rms errors for the case of advanced detectors are 10 times less than for the initial detectors. 
This is a result of simplified assumptions for the spectral density of shot noise for the advanced detectors, 
namely it is taken exactly 100 times less than for the initial detectors (for frequency range of 500 Hz to 
1 kHz the shot noise dominates). 

The GEO600 detector in narrowband configuration tuned to a certain frequency (1 kHz in this case) 
is able to achieve accuracies of estimation of parameters in the phase of the signal (position on the sky, 
frequency and spindown parameters) several times better than kilometer-length initial detectors. This 
is done at the expense of loss of information at other frequency bands. We see that the accuracy of 
the estimation of the amplitude h Q and the angle 9 for the narrowband case is less than for any other 
detector. This is because the narrowband GEO600 detector will see mainly one component of the signal 
{h,2 in this case) where parameters h and 9 degenerate into one parameter. The information about the 
other component is small resulting in a small accuracy for joint estimation of parameters h a and 9. It 
implies that for the narrowband GEO600 detector it is more suitable to model the signal as consisting 
only of the second Q12) component. Then the parameters h a and 9 merge into an effective amplitude 
parameter h' a — h Q sin 2 9 [cf. Eqs. (§])-(§])]. We have performed the Monte Carlo simulations for this case. 
On the left top panel in Figure 1 we have shown the cumulative distribution function for the relative rms 
error of the amplitude h' a . We have checked that the position error Af2 and the spindown relative errors 

<y(fo)/\fo\ (k = 0, ... ,4) are almost the same as for the case of the two-component signal model indicating 
negligible contribution of the first component [hi) to the information about the phase parameters. 

From Table 2 we see that the initial laser-interferometric gravitational-wave detectors will be able to 
locate the strongest fast and young neutron star gravitational-wave emitters at a distance of 1 kpc with 
accuracy of order 10~ 7 sr and estimate the relative accuracies of frequency (k = 0) and four spindowns 
(k = 1, ... ,4) within a factor of order 10~ 2 ( 4 ~ fe '. Amplitude can be estimated within a few percent and 
angle 9 within a few percent of a radian. 

In Figures 2 and 3 we have performed the Monte Carlo simulations of the rms parameter errors for the 
eight models of the phase summarized in Table 1 above. Figure 2 is for long observation time T, = 120 
days and Figure 3 is for short observation time T D — 7 days. Because the amplitude of the signal for our 
model increases as square of the frequency [see Eq. (||)] the rms errors of the parameters for fast rotating 
stars are smaller than for slowly rotating ones. For the case of young neutron stars the spindowns are 



larger than for the old ones [see Eq. (|15|)] and this leads to more spindown terms in the phase for young 
stars in order to meet the 1/4 of a cycle criterion (see Table 1). Since spindown parameters are strongly 
correlated (see Appendix C) increasing their number decreases the accuracy of their estimation. Thus 
accuracies for estimation of parameters for young neutron stars tend to be less than for the old ones. 

(4) 

We observe that the accuracy of estimation of the phase parameters (except for the fourth spindown f a ) 
is in all cases very good. The distributions of the rms errors of the amplitude parameters h a and 9 for 
various distinct models of the phase in the case when the frequency f and consequently the amplitude 
h a is the same are indistinguishable (see top panels of Figures 2 and 3). This shows that the amplitude 
parameters and the phase parameters in our signal model are effectively uncorrelated. 

In Figures 4-6 and Table 3 we have presented results of the same simulations as in Figures 1-3 and 
Table 2 but for directed searches. Here because there are two unknown parameters less in the phase of the 
signal (the right ascension a and the declination S are assumed to be known) the accuracy of estimation 
of the phase parameters increases. For the initial detectors the relative rms errors of the frequency 
(k = 0) and the first two spindowns (k = 1, 2) are of the order of 10~ 3 ( 4 ~ k \ and the relative rms errors 
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7.9xl0~ 2 
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38 
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Table 2: The quartiles (To) for the Monte Carlo simulated distributions of the rms errors of the signal 
parameters for the individual detectors in the case of all-sky searches. The observation time T a = 120 
days. The neutron star parameters are the same as in Figure 1, the frequency f a = 500 Hz, and the 
spindown age t = 40 years. The model of the signal's phase is described by si = 4, S2 = 3, and S3 = (cf. 
Table 1). The dimensionless amplitude of the waveform is h Q = 1.1 x 10~ 23 . For the GEO600 detector 
we use two noise curves: wideband and narrowband tuned to 1 kHz with the bandwidth of 30 Hz. The 
last column shows the quartiles of the Monte Carlo simulated distributions of the signal-to-noise ratio d 
[given by Eq. (p8|)], it is taken from Table II of Paper I. 

of the third and the fourth spindown are of the order 10 -4 and 10 _1 , respectively. This gives a decrease 
with respect to the all-sky searches by factors of 10 4 , 10 3 , 10 2 , 10 2 , and 10 for the frequency, the first, 
the second, the third, and the fourth spindown respectively. Again the distributions of the rms errors 
of the amplitude parameters h a and 9 for the models with the same frequency f a are indistinguishable 
(see top panels of Figures 5 and 6). Moreover the distributions in the top panels of Figures 2 and 3 
(all-sky searches) are indistinguishable from the corresponding distributions in Figures 5 and 6 (directed 
searches). Since the signal model for all-sky searches and the corresponding model for directed searches 
differ only in the number of the unknown parameters in the phase of the signal this confirms an effective 
decorrelation of the amplitude and the phase parameters. This fact will enable us (in Sec. 6 below) to 
give an assessment of the rms errors of the phase parameters independently of the physical mechanisms 
generating the gravitational waves from a spinning neutron star. 

4 Estimation of the proper motion of the neutron star 

To describe the proper motion of the neutron star we shall introduce the parameters that are customarily 
used by the pulsar astronomers. When the star moves with respect to the SSB reference frame its right 
ascension a and declination S change. To first order we can assume that they change linearly with time, 
i.e. we have 

a = a + fi a t, 5 = S + figt (30) 

Let us denote by v ns the velocity of the star w.r.t. the SSB reference frame and by v ns j_ the component 
of the velocity v ns perpendicular to the unit vector no in the direction of the star (no is along the line of 
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Table 3: The quartiles (To) for the Monte Carlo simulated distributions of the rms errors of the signal 
parameters for the individual detectors in the case when the position of the source in the sky is known. The 
observation time T = 120 days. The neutron star parameters are the same as in Figure 1, the frequency 
jo = 500 Hz, and the spindown age r = 40 years. The model of the signal's phase is described by s\ = 4, 
s-2 = 3, and S3 = (cf. Table 1). The dimensionless amplitude of the waveform is h a = 1.1 x 10 -23 . 
For the GEO600 detector we use two noise curves: wideband and narrowband tuned to 1 kHz with the 
bandwidth of 30 Hz. The last column shows the quartiles of the Monte Carlo simulated distributions of 
the signal-to-noise ratio d [given by Eq. (f28|)], it is taken from Table II of Paper I. 



sight to the star at t = for the observer located in the SSB). To first order in time it is easy to express 
the velocity v ns ^ in terms of the parameters \x a and /ig. We have 

— n a sinao cos Sq — \ig cosao sin<5o \ 

fi a cos ao cos Sq — fig sin ao sin Sq . (31) 
Us cos S J 

The total proper motion of the star fj, := |v ns j_|/r is given in terms of fj, a and fig as follows: 



/i = yJ/4 + til COS^o. (32) 

In our Monte Carlo simulation we have considered the extreme case of a nearby neutron star at a 
distance of 40 pc. We have taken the transverse velocity |v ns jj of the star to be 10 3 km/s. Careful 
modelling of the pulsar statistics has shown that the pulsar velocities extend to these values |11, L|. For 
|v n sj_| = 10 3 km/s and r$ = 40 pc the contribution ( Jilf ) to the phase due to the proper motion is ^4 
cycles for frequency f = 500 Hz (this value we have used in the simulation) and it decreases linearly 
with distance. 

One can show that in the extreme case we consider here the phase model consistent with the 1/4 of a 
cycle criterion for 120 days of observation time and young neutron star with spindown age r = 40 years 
reads [cf. Eqs. @ and @] 



•w - + *±T.£m + f »» ■ - w tK + 1 (« • + i /.. (33) 
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For the phase given by Eq. ( |33| ) the two-component gravitational-wave signal defined by Eqs. (|l|)-(^) 

(1) (4) 

depends on the following 14 parameters: h Q , ao,/j, a , 8q, Us, i>, 0, $0j fo, fo, ■ ■ ■ , fo- 

In the Monte Carlo simulation for each set of randomly generated angles {ao, So, tp, l, 8} we have also 
randomly chosen the direction of the star's motion on the sky keeping the total proper motion /i of the 
star fixed and equal to (10 3 km/s)/(40 pc) = 5.3 x 10 3 mas/yr (mas/yr stands for milliarcseconds per 
year, units usually used in pulsar astronomy). 

The results of the numerical simulations are summarized in Figure 7. We have assumed Hanford LIGO 
detector in the advanced configuration and 120 days of observation time. The extreme case considered 
here gives interesting accuracies of the estimation of the proper motion parameters resulting in the median 
of the distributions of a{fjL a ) and ct(iis) to be ~ 2 x 10 3 mas/yr (i.e. 40% rms error). As the distance to 

(4) 

the source is small the estimation of the fourth spindown parameter f a with a good accuracy of ~10% is 
also possible. 

5 Simplified models of the gravitational-wave signal 

The phase of the gravitational-wave signal from a spinning neutron star changes at the detector on the 
Earth over a characteristic time of less than 100 ms whereas the amplitude of the signal changes over one 
day. This means that the detection of a long continuous signal requires an accurate model for its phase. 
As shown in Appendix A even 1/4 of a cycle difference between the filter and the signal can lead to loss 
of the signal-to-noise ratio by 10%. 

Accurate modelling of the amplitude is not so crucial. Consequently the information on the phase 
parameters in the amplitudes of the signal is much smaller than in the phase. This can clearly be seen 
from our simulations presented in Sec. 3 which show that the phase parameters can be estimated much 
more accurately than the amplitude parameters (see Figures 1 to 6). As a result the phase parameters 
effectively decorrelate from the amplitude parameters. This decorrelation is shown in our numerical 
simulations (see discussion at the end of Sec. 3 and Figures 2, 3, 5, 6, see also Figure 15). Therefore we 
can expect to obtain a model of a signal with constant amplitudes (i.e. with no amplitude modulation) 
that reproduces the covariance matrix for the phase parameters of the full signal. Such a model is 
presented in Sec. 5.1 below. 

We also find it useful to have a linear parametrization of the phase of the signal i.e. a set of parameters 
such that the phase is a linear function of the parameters. One advantage of the linear parametrization 
is that the maximum likelihood estimators of the parameters are unbiased |0. The other property of a 
linearly parametrized phase is that the covariance matrix is independent of the values of the parameters. 
In real data analysis schemes optimization of the codes when using a linear parametrization would be 
independent of the values of the parameters, i.e. of the parameter region searched. We find that the 
linear parametrization is not possible in general. In Sec. 5.2 we show however that it can be achieved 
approximately by neglecting a term in the phase arising from Earth diurnal motion and all contributions 
from the spindowns in the Earth spin and orbital terms. 

In Sec. 5.3 we present another linear model (called linear model II) that has been investigated by 
one of us . It is obtained by neglecting completely Earth diurnal motion and all spindowns in Earth 
orbital motion terms. The simplest, polynomial in time phase model, is introduced in Sec. 5.4. In this 
model the motion of the detector w.r.t. the SSB is ignored. In Sec. 5.5 we describe the results of Monte 
Carlo simulations of the covariance matrix to compare the simplified models of the signal with the exact 
two-component model given in Sec. 2. 

5.1 Constant amplitude model 

The constant amplitude model of the gravitational- wave signal has the phase of the exact model whereas 
the time dependent amplitudes are replaced by some constant effective values. We have obtained these 
effective values from the analytic formula for the signal-to-noise ratio given in Appendix B of Paper I. 
The signal-to-noise ratio consists of a term proportional to the square root of the observation time and 
some oscillatory term. For the observation times longer than several days the term proportional to the 
square root of the observation time strongly dominates. Our effective amplitudes take into account only 
this dominant part of the signal-to-noise ratio. Thus for the two-component signal h that we consider 
the constant amplitude model depends on the two constant amplitudes h \ and h i as well as on the two 
initial phase parameters $oi an d $02'- 

h(t) = hoi sin + *qi] + h o2 sin [2*(t) + $02] , (34) 
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where the phase $ coincides with the exact phase model given by Eq. ( |12| ) . The constant amplitudes h a \ 
and h o2 read [cf. Eqs. (85), (86), and Appendix B in Paper I] 



where 



h i = h a sin C| sm29\y / F 1 (L)e 1 (5) cos 4-0 + Gi(t)e 2 (6), (35) 
h o2 = h sin ( sin 2 9^/F 2 {i)e 1 {5) cos 4^ + G 2 (i)e 2 (5), (36) 

Fi(t) = -^sin 4 t, F 2 (l) = isin 4 /,, 

Gift) = — sin 2 t (l + cos 2 t) , G 2 (t) = - (l + 6 cos 2 t + cos 4 t) , 
16 4 

ei(S) = 4ji cos 4 (5, e 2 (<5) = 4j 2 — 33 cos 25 + ji cos 2 25, 

h = — (4 -20 cos 2 A + 35 sin 2 2 7 cos 4 A) , 
256 

j 2 = (68 - 20 cos 2 A - 13 sin 2 2 7 cos 4 A) , 

j a = (28 - 44 cos 2 A + 5 sin 2 2 7 cos 4 A) . 

128 

5.2 Linear model I 

The signal h! of the linear model I has two constant amplitudes given by Eqs. ([55]) and ( |36| ) so it can be 
written as 

h'(t) = h i sin [$'(0 + $ i] + h 2 sin [2$'(t) + $02] • (37) 

The phase $' of the signal /i' differs from the exact phase model $ given by Eq. ( |l2| ) by neglecting 
all spindowns in the phase modulation due to the orbital motion of the Earth and by discarding the 
component of the vector defining the position of the detector w.r.t. the SSB which is perpendicular to 
the ecliptic (see Sec. 2. 2 of Paper I). The function <E>' is thus given by 

* ( fc ) fk+l 

= 2vr]T/ 



fe =0 + 
2nfo 

H { (cos e sin a cos 5 + sine sin 8) [Res sin (6 + f2 D t) + Re cos A cose sin (6 r + Q r t)] 

c 

+ cos a cos S [Res cos (<fi + &ot) + Re cos A cos (0 r + O r t)]} . (38) 
After introducing the two new parameters 

a i '■= fo (cos e sin a cos 5 + sine sin 5) , a 2 := f a cos a cos 6, (39) 

r— . (i) (•) 

the phase ( p8[ ) becomes a linear function of the new parameters f , f Q , . . . , f , a,\, and a 2 . 

5.3 Linear model II 

For the linear model II the signal h" is defined by 

h"(t) = hoi sin [$"(*) + $01] + h o2 sin [2$"(t) + $02] , (40) 



where again the constant amplitudes h Q \ and h o2 are given by Eqs. (35) and (p6|). In the phase of 
the signal compared with the exact phase model $ given by Eq. (|l2) we neglect all spindowns in the 
phase modulation due to the orbital motion of the Earth and we discard the whole modulation due to 
the Earth's diurnal motion. The function <E>" is thus given by 

s (fc) j.fe+1 

*"«> - 

27r/ 

H -Res [cos a cos 6 cos {<j> + fl a t) + (cos e sin a cos S + sine sin S) sin [6 + r2 £)l . (41) 

c 

As in the previous model after introducing the parameters 

a i '■= fo (cos e sin a cos (5 + sin e sin (5) , a 2 := / Q cos a cos (5, (42) 

. — . _ . <i) («) 

the phase ( [41[ ) becomes a linear function of the parameters f a , fo: ■ ■ ■ > fo, c^ij and ol 2 . 

A version of the above model including only the second component was used in a simplified study of 
data analysis of gravitational- wave signals from spinning neutron stars made by one of the authors |lE|] . 

Let us observe that both the constant amplitude model of Sec. 5.1 and the linear model I of Sec. 5.2 
coincide with the linear model II for a detector located at the north or the south pole (then A = ±90°). 
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5.4 Polynomial model 

In this model the signal h'" is given by 

h'"(t) = h ol sin + $ i] + h o2 sin [2$'"(i) + $ 02 ] , (43) 

with the constant amplitudes and /i e 2 defined by Eqs. (^) and (|36|). In the phase <&'" of the signal 
( [43| ) we discard all the terms due to the motion of the detector relative to the SSB. This leads to the 
following polynomial in time model of the phase 

* (fc) fk+l 

= ^E/oprnjr (44) 

fe=0 v ' 

For the case of directed searches i.e. when the position of the source in the sky is known, the phase 
is a linear function of the unknown parameters: initial phase, frequency, and spindowns and thus the 
exact phase model is linear. In Sec. 5.5 below we find that for directed searches the polynomial model 
reproduces accurately the covariance matrix for the exact model. We investigate the polynomial phase 
model in more detail in Appendix C. 



5.5 Simplified models vs. the exact model 

We have carried out the Monte Carlo simulations to compare our simplified signal models with the exact 
model given in Sec. 2. In our comparisons we have considered both all-sky searches and directed searches 
and we have chosen two observation times: 7 days and 120 days. We have found that the constant 
amplitude model of Sec. 5.1 reproduces very accurately the results of the exact model (see Figures ^| 
and ^| for results of Monte Carlo simulations) . Also the linear model I reproduces the rms errors of the 
exact model well (see Figures || and |^). We have compared the percentiles of the distributions of the 
rms errors of the phase parameters for this model with the percentiles of the distributions of the exact 
model. We have found that for all-sky searches and observation time of 7 days, frequency f Q = 500 Hz, 
and spindown age r = 40 years, the maximum differences for various parameters range from ~15% to 
~60%. For observation time of 120 days and the same frequency and spindown age the corresponding 
differences range from ~7% to ~50%. For the linear model II the deviations of the rms errors from the 
exact model are large (see Figures || and ^|) . 

We have found that for directed searches the polynomial phase model introduced in Sec. 5.4 reproduces 
the results of the exact model very well. For observation time of 7 days the maximum percentile differences 
of the distributions of the phase parameter errors are around 5% and for observation time of 120 days 
they are less then 1%. 

Numerical calculations of the covariance matrices for the constant amplitude models are enormously 
simplified compared to the amplitude modulated signals. Since the model of the phase is independent of 
the physical mechanisms generating the gravitational radiation, the simpler constant amplitude models 
characterize well a general continuous gravitational-wave signal and they can be a useful tool in theoretical 
analysis of the gravitational-wave signals from spinning neutron stars. However we cannot compromise 
the detectability of signals for simplicity of the constant amplitude and linear phase models and we do 
not insist on using them in real data analysis schemes. Nonetheless for searches over certain limited 
parameter space or for shorter observation times the simplified models (both constant amplitude and 
linear phase models) may be useful. 



6 Dependence of the covariance matrix on the observational 
parameters 

It is clear that covariance matrices of various signals considered in the present paper depend on the 
observation time T a . We also notice that the phases of our signals and consequently the covariance 
matrices depend also on the deterministic phases <j> , 4> r and on the latitude A of the detector's location 
[see Eqs. (p^)-(p^)]. The values of the phases <p and <j> r depend on the initial time of the observation. It 
is important to know these dependencies and therefore we study them in detail in this section. 

The results of the previous section indicate that the constant amplitude linear phase model of Sec. 
5.2 reproduces well the covariance matrices of the exact model. For directed searches even the simpler 
polynomial model of Sec. 5.4 is adequate. Therefore we adopt these two models to study the dependence 
of covariance matrices on T a , <p , </> r , and A. Moreover we restrict ourselves to the one-component versions 
of these models. 
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Both models of Sec. 5.2 and of Sec. 5.4 are the constant amplitude models. We first derive the general 
formulae for the covariance matrix valid for the constant amplitude signal with any phase model, provided 
the signal is narrowband around some frequency f a . Such a one-component signal can be written as 

h(t;h ,C) = /i sin*(t;C), (45) 

where h is a constant amplitude and £ denotes parameters entering the phase 'J of the signal. We collect 
all the signal parameters into the vector — (h a , £). Assuming that over the bandwidth of the signal the 
spectral density Sh of the noise is nearly constant and equal to Sh(fo) and that during the observation 
time the phase \t has many cycles the optimal signal-to-noise ratio d :— ■J (h\h) for the signal ( f45| ) equals 

, ^ h \fT 



(46) 



The Fisher matrix T for that signal can be written as [cf. Eq. ([19)]: 

2 [ T ° /2 dh dh 

Teie * = wJ-T o/ 2dOido; dt ' (7) 

where T a is the observation time and the observation interval is [— T /2, T /2\. The components of the 
matrix T read: 

(48) 
(49) 

CiO > (50) 
where the reduced Fisher matrix T is defined by 

1 [ T ° /2 9* 5* 





d 2 


" h? 




S 0, 




d 2 f, 



T j-T a /2 d(i d(j ^ ^ 

The covariance matrix C approximated by the inverse of the Fisher matrix T has the components 

(52) 
(53) 
(54) 





h 2 


Ch h 




Chad 


s o, 















CiCj 



<fc) 

It is sometimes convenient to replace the spindowns f a by the dimcnsionlcss parameters ujk defined as 

- : = JkTTy. f ° T ° +1 > ^ 

ft) 

where T is the observation time. The rms errors of the spindown parameters f a are related to the rms 
errors of the dimensionless parameters u>k by 



(*o (k + iy i 

°{fo)= 1 I ' =fc+r ^K , fc = 0,..., S . (56) 

27T To 



The significance of the parameters can be explained as follows. In the polynomial phase model 
(discussed in Appendix C) the covariance matrix for the parameters Uk is independent of the observation 
time T Dl so it is completely determined by the number of spindowns included in the phase. In all more 
complicated phase models considered in this paper, including the exact phase model introduced in Sec. 2, 
the phase is the sum of the polynomial in time part and extra terms due to the motion of the detector w.r.t. 
the SSB. The polynomial part gives dominant contribution to the number of cycles in the phase and it 
roughly determines how the covariance matrices depend on the observation time T Q . Using the parameters 

u>k instead of f a we absorb the polynomial phase model dependence on T Q into the very definition of the 
parameters ujk- If the covariance matrix for the parameters u>k in a non-polynomial phase model depends 
on T , this dependence is a measure of how this model is different from the polynomial phase model. 

For all-sky searches we have studied the dependence of covariance matrices on T , <f> r , <fi , and A using 
the constant amplitude linear phase signal of Sec. 5.2. We have rewritten the phase of this signal in terms 
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of the parameters ( J38| ) and (|55|) (then the phase is a linear function of these parameters) . The signal can 
be written as [cf. Eq. @] 



h(t-KX) = /i sin*(t;C), (57) 
s f t \ k+1 2-k 

*(t;C) = $o + ^ a)fc I — 1 H {ai [Res sin (0 O + o t) + Re cos A cos e sin (</v + ri r t)] 



■.To J c 

k=0 v 7 

+a 2 [-Res cos (0 o + Q a t) + Re cos A cos (4> r + £l r t)]} , (58) 

where C — (ho, a.%, u>o, . . . , co s ). Thus the signal ( |57j) depends on 5 + s parameters. We have 

computed the covariance matrices for this signal by means of Eqs. (|5^)-(|54|). The results are discussed 
below. 

To study the dependence of the rms errors of the phase parameters on the observation time T D we 
have considered a wide range of observation times from 1 hour to 4000 days. We have taken the latitude 
A of the LIGO detector in Hanford and we have assumed the signal-to-noise ratio d — 10. We have 
considered five phase models given by Eq. (p8|) with s = 0, . . . , 4 spindowns. The results of the numerical 



calculations are shown in Figure 10. We see that for observation times T up to ^20 days the rms 
parameter errors decrease. For observation times from ^20 days to ~100 days the rms errors for initial 
phase $o and parameters ol\ and a 2 tend to decrease. The rate of the decrease is the smaller the more 
spindown parameters enter the phase of the signal and for higher spindown models the errors stay even 
constant for certain ranges of the observation time, whereas the frequency and the spindown errors in 
these ranges tend to increase. For high spindowns there is a range of observation time for which rms errors 
stay constant. For observation times from ~100 days to ~1000 days all the rms errors decrease and for 
observation times greater than ^1000 days the rms errors level out to constant values. We have verified 
that all these constant values are to a very good accuracy equal to the rms errors of the polynomial 
phase model for the corresponding number of spindown parameters (these constant values can be found 
in Table 4 of Appendix C). Thus for observation times of more than ^3 years the effect of the motion of 
the detector relative to the SSB has negligible effect for the rms errors of the initial phase, the frequency, 
and the spindown parameters. 

The dependence of the rms phase parameters errors on the phase 4> r which is determined by the 
initial position of the Earth in its diurnal motion is small. The errors vary by no more than ^10% for 
the observation time T Q = 7 days and for the phase models (^8|) with s = 0, . . . , 4 spindowns (see Figure 
O, where the model with s = 4 spindowns is studied) assuming the signal-to- noise ratio d = 10. This 
dependence weakens when the observation time T Q increases. 

The variations of the rms errors of the phase parameters with the initial phase 4> a of the Earth's 
orbital motion are larger and depend in a complicated way on the observation time T Q and the number 
s of spindowns included in the phase model. For the observation time T a — 120 days and the phase 
model (|5^) with s = 4 spindowns (assuming the signal-to-noise ratio d = 10) the variations are less than 
^10% for all parameters except for parameters a\ and a 2 where they are of order 100% (see Figure |l2| ). 
For the observation time T Q — 120 days and the signal-to-noise ratio d = 10 the phase models ( |58| ) with 
s = 0, . . . , 3 have the rms errors of all parameters that vary less with (f> compared to the s = 4 model, 
with the exception of parameters u\ and a 2 , for which the variations range from ~500% to ^900%. For 
the observation time T a — 7 days and the phase models (|58|) with s = 0, . . . , 4 spindowns (assuming the 
signal-to-noise ratio d = 10) the variations are less than ^10% for all the parameters except again for the 
parameters a± and a 2 . For these parameters the variations range from ~10~ 7 % for s = 4 up to ^10 3 % 
for s = 0. 

The rms errors of the phase parameters exhibit a sharp increase when the detector's location ap- 
proaches the south or north pole. The dependence is weaker for the phase models ([38]) with smaller 
number s of spindowns included and it gets stronger when the observation time T a increases. However 
the latitudes of the laser-interferometric detectors currently under construction range from 30.56 to 52.25 
degrees (see Table I in Paper I) and in this interval the rms errors vary by no more than ^30% for the 
observation time T D = 120 days and the phase model (H|) with s — 4 spindowns (see Figure |l3|) assuming 
the signal-to-noise ratio d — 10. The sharp increase of the rms errors when A tends to ±90° explains why 
the linear phase model II of Sec. 5.3 gives the errors significantly larger compared to another constant 
amplitude models (see Figures 8 and 9). 

In Table 4 we have collected the rms errors of the initial phase <J>0j the parameters a±, a 2 , and the 
(fc) . — . 
spindowns f Q [k = 0, . . . , 4) in the case of all-sky searches. We have considered five phase models (|38|) with 

s = 0, . . . , 4 spindowns included and the observation times T Q of 7 and 120 days. The errors are calculated 

for the signal-to- noise ratio d = 10 and they are inversely proportional to d. They are independent of the 

values of the spindown parameters in the signal. As a reference we quote the values of the spindowns for 

the Crab pulsar: f a = -3.773 x lO" 10 s~ 2 , f = 0.976 x 10" 20 s^ 3 , and f a = -0.615 x 10~ 30 s~ 4 . We 
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see that in many cases the value of the rms error of a certain parameter does not change if we go from 
a model with k spindowns to a model with k + 1 spindowns. For example for observation time of 7 days 
the rms error of frequency is almost the same for all the models. This merging of the errors is visible 
clearly in Figure 10 and occurs for certain ranges of the observation time. It is related to the merging for 
the polynomial phase model studied in Appendix C (see Figures 14 and 15). 

The accuracy Af2 of the position of the source in the sky can be expressed in terms of the rms errors 
of a±, a.2, and f a by means of the rule of propagation of errors. 

We found in the previous section that for directed searches the polynomial phase model is adequate. 
Consequently in the case of directed searches we can completely neglect the dependence of the covariance 
matrices on the initial phases <p r and <j) and also on the latitude A of the detector's location. The rms 
errors of the initial phase $o and the dimensionless spindown parameters cut (k = 0, . . . , 4) can be found 
in Table 4 of Appendix C. They are constant numbers, independent on the observation time T Q . However 

(JO 

the rms errors of the spindown parameters f a (k = 0, . . . , 4) do depend on the observation time T Q and 
this dependence is described by Eq. ([56]). 

(*) 

In Table 5 we have given the rms errors of the initial phase $o and the spindowns f a (k = 0, . . . , 4) 
for directed searches. We have considered five polynomial phase models with s = 0, ... ,4 spindowns 
included and the observation times T Q of 7 and 120 days. The errors are calculated for the signal-to-noise 
ratio d — 10 and they are inversely proportional to d. 

The fact that in some cases (both in all-sky and directed searches) the errors of the parameters do 
not change when we go to a model with one more spindown is a result of a special choice of the initial 
time of the observation which was chosen to be in the middle of the observation interval. This effect is 
studied in detail in Appendix C (see Figure 14). 
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Table 4: The rms errors of the initial phase $o, the parameters a\, and the spindowns f Q (k = 0, . . . , 4) 
for all-sky searches. We have approximated the gravitational-wave signal by the one-component constant 
amplitude and linear phase model of Sec. 5.2. We have considered five phase models with s = 0, ... ,4 
spindowns included and the observation times T a of 7 and 120 days. The signal-to-noise ratio d — 10. We 
have assumed the latitude A of the Hanford LIGO detector, we have also put <j> r — 1.456 and <j) = 0.123. 
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Table 5: The rms errors of the initial phase <&o and the spindowns f (k = 0, . . . , 4) for directed searches. 
We have approximated the gravitational-wave signal by the one-component constant amplitude and poly- 
nomial phase model of Sec. 5.4. We have considered five phase models with s — 0, ...,4 spindowns 
included and the observation times T a of 7 and 120 days. We have assumed the signal-to-noise ratio 
d= 10. 

A 1/4 of a cycle criterion 

The fitting factor FF between a signal h — h(t; 6) and a filter h! — h'(t; 6') (6 and 6' are the parameters 
of the signal and the filter, respectively) is defined as |13| 

FF := max ' (59) 

e ' y /(h(t;0)\h(t;O))yJ(h'(t;6')\h'(t;0')) 

1 — FF gives the fraction of the signal-to-noise lost when using a filter not perfectly matched to a signal. 
For narrowband signals around the frequency f the scalar product (-|-) can be computed from the formula 



H'y\l>-.) =5 [ ° hx(t)h»(t)dk, (60) 

-T /2 



Sh(fo) 



where Sh is the one-sided noise spectral density and T a is the observation time. 
Let us assume that the signal and the filter can be written as 

h(t;8) = hoSm^(t;C), (61) 
h'(t;d') = h sm*'(t;C'), (62) 

where h and h' a are constant amplitudes, £ and C' denote the parameters entering the phases ^ and 
of the signal and the filter, respectively. We substitute Eqs. ( |6l| ) and ( p2|) into Eq. (|59|). Using Eq. ( |60|) 
we obtain 

i r T ° /2 

FF S max — / cos C) - C )] dt. (63) 

C T-o J-T /2 

The fitting factor attains its maximum value of 1 when the functions $ and 'J'' are the same and when 
values of the parameters in the signal and the filter coincide. When "J and W are not the same because in 
the phase of the filter we have not taken into account some effects present in the signal the fitting factor 
is less than 1. Moreover the values of the parameter in the filter that maximize the correlation integral 
will be biased, i.e. shifted away from the true values in the signal. 

For the simplest nontrivial case of the difference between the phase of the filter and the phase of 
the signal consisting of a constant term (constant phase <f) and a term linear in time (with a constant 
frequency /) the fitting factor (B3h equals 



1 



T /2 



FF S max — / cos (2nft + 0) dt. (64) 

<t> T J^To/2 



We easily get 

11= max 



COS( 



(65) 
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If the frequency / is such that it produces not more than 1/4 of a cycle during the observation time T 
(i.e. fT Q < 1/4) then from Eq. @ we obtain 

- sin ? 
FF > = 0.900. 

4 

Thus characterization of the 1/4 of a cycle criterion is the following: for the simplest nontrivial case 
discarding a term in the phase that contributes less than 1/4 of a cycle over observation time leads to a 
loss in signal-to-noise ratio of not more than 10%. 



B Number of cycles 

The model of the phase of the signal at the detector given by Eq. ( |l2| ) can be rewritten as 

Sl S 2 S 3 



(66) 



fc=0 



k=0 



k=0 



(o) (r) 

where Nk, Nf, and denote the numbers of cycles arising from the term polynomial in time, the 
orbital motion and the rotational motion of the detector, respectively: 
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(67) 



(!=) 



where f g and f g are respectively the frequency and the spindown parameters of the gravitational-wave 
signal. The moduli of the quantites Nk, N^ and N^ can be estimated as follows: 
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(68) 
(69) 



We estimate the maximum values of the spindown parameters f g from the following formula 



(70) 



where r is the spindown age of the neutron star. Using Eq. ( fTOj) the quantities N k , Nj^ and N^ r> can 
be estimated as follows: 
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(71) 
(72) 



Let us note that the ratios of the different contributions for the fixed k do not depend on the gravitational- 
wave frequency f g and the spindown age r: 
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(73) 



From the formulae given above it is easy to obtain how many terms are needed in each of the three 
series on the right-hand side of Eq. ( |l2| ) for the eight models of the phase considered in our simulations 
in order to meet the criterion that all terms that contribute more than 1/4 of a cycle are included. The 
results are summarized in Table 1 of Sec. 2. 



C Polynomial phase model 

The ability to control the values of the Fisher information matrix for the frequency of the signal and 
its derivatives may be useful in real data processing schemes. Therefore in this appendix we study the 
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dependence of the Fisher matrix on the choice of the instant of time at which the instantaneous frequency 
and spindown parameters are defined. 

We use the one-component signal with a constant amplitude and the phase polynomial in time, which 
is introduced in Sec. 5.4. Such a signal can be written as 



h(t;h a ,C) = /i sin*(i;C), 



(74) 



where h a is a constant amplitude and £ denotes the parameters entering the phase \& of the signal. As 
the signal (Q) has the constant amplitude, the formulae (|^|)-([54|) from Sec. 6 can be applied here. But 
here unlike in Sec. 6 we choose the observation interval to be [0,T o ], where T a is the observation time. 
We also introduce the dimensionless variable x := t/T a . Then the formulae ([52|)-([54]) for the components 
of the covariance matrix C can be rewritten as 

h 2 

(75) 



C t .t< 



d? 1 

0, 



(76) 
(77) 



where the optimal signal-to-noise ratio d is given by Eq. ([16|) and the reduced Fisher matrix T is defined 
by 



■ dx. 



Using the dimensionless variable x the polynomial phase is defined as 



(78) 



*0r;C) = $o+ }u k (x-x ) 



k+l 



fe=0 



(79) 
, s) are the 



where xq is an arbitrarily chosen initial time, $0 denotes the initial phase, and ujk (k = 0, 
spindown parameters, so £ = ($0, ^o, . . . , u> s )- 

The reduced Fisher matrix T for the polynomial phase model ( |79| ) has the structure very similar to 
that of the (s + 2)-dimensional Hilbert matrix. For xq = T is the (s + 2)-dimcnsional Hilbcrt matrix. 
For s — 4 (four spindown parameters included) the reduced Fisher matrix T computed for an arbitrary 
xq equals 
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where 
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After putting xq = into the above formulae one obtains the 6-dimensional Hilbert matrix: 
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Table 6: The rms errors <r(Ck) of the phase parameters C = (^o, ^o, ■ ■ • ,1^4) for the polynomial models 
with s = 0, . . . , 4 spindowns. We have assumed the initial time xq = 0.5. 



We have computed the rms errors of the phase parameters for five different phase models (|7^) with 
s = 0, . . . , 4. The rms errors cr(Cfc) are the square roots of the diagonal elements of the reduced covariance 
matrices C = T^ 1 . In Figure |l4| the errors cr((k) are plotted as functions of the initial time xq. It is 
seen from the figure that the errors depend on the choice of the initial time xq. This can be explained 
as follows. The phase ^ of the signal at a given time x does not depend on the choice of the initial time 
xq ■ Thus for the two different choices xqi and X02 of the initial time xq we have 
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(82) 



where the parameters $oi> w fei and $02, Uk2 are associated with the initial times xoi and Xq2, respectively. 
Using Eq. ( p2[ ) one can, for fixed s, express the parameters £i = (<&oi> ^oii • • • j w si) as linear combinations 
01 C2 = (*02, ^02, ■ ■ • ,w s2 ): 
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where are the components of the transformation matrix J given by 
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(84) 



Thus the change of the initial time from xqi to X02 corresponds to estimating a different set of parameters 
given by a linear combination of the original parameters. For example the new frequency parameter 
defined as the first derivative of the phase at time X02 (divided by 2ir) is a linear combination of the 
frequency parameter and spindown parameters at time Xo%. Let us denote by C(xqi) and C(xq2) the 
covariance matrices for the parameters Ci; and C21, respectively. It is easy to see that these matrices are 
related by 



C(x m ) = JC(x 01 )J J 



(85) 



The dependence of the rms errors of the signal parameters on the choice of the initial time for the 
full phase model is very similar to that observed in the polynomial model considered above. This is 
illustrated in Figure |l5|, where we have plotted the simulated cumulative distribution functions of the 
parameter errors for the full model of the signal described in Sec. 2. The initial time t a — used in the 
simulations coincides with the middle of the observation interval which is [— T G /2, T a /2], so it corresponds 
to xo — 0.5. The comparison of Figures 14 and [15] shows that the curve merging observed in Figure [l5| 
exactly corresponds to intersection points observed in Figure [l4]for xo — 0.5. 

Finally in Table § we give the rms errors (diagonal components of the inverse of the Fisher matrix) 
of the parameters C m various polynomial models. We consider models with s = 0, . . . , 4 spindowns and 
we assume xq — 0.5. 
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Figure 1: Cumulative distribution functions of the simulated rms errors of the signal parameters for 
the individual detectors in the case of all-sky searches. The observation time T — 120 days. We 
assume that star's ellipticity e = 10~ 5 , its moment of inertia w.r.t. the rotation axis I = 10 45 g cm 2 , its 
distance from the Earth r Q = 1 kpc, and the frequency f = 500 Hz [these values give the dimensionless 
amplitude of the waveform h Q = 1.1 x 10~ 23 , cf. Eq. (|B|)]. The spindown age r = 40 years. The model 
of the signal's phase is described by s\ — 4, S2 — 3, and S3 = (cf. Table 1). The lines given on the 
diagrams correspond to various detectors: advanced Hanford (dotted/double dashed), initial Hanford 
(solid), VIRGO (dotted), wideband GEO600 (dashed), narrowband GEO600 (double dotted/dashed), 
and TAMA300 (dotted/dashed). On the left top panel we have also put the extra double dotted/double 
dashed curve. It describes cumulative distribution function of the relative rms error of the amplitude 
h'o = h sin 2 9 of the second component of the signal as measured by the GEO600 narrowband detector. 
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Figure 2: Cumulative distribution functions of the simulated rms errors of the signal parameters for 
various models of the signal in the case of all-sky searches. The observation time T a = 120 days. The 
initial Hanford detector is assumed in the simulations. The neutron star parameters are the same as in 
Figure 1. The lines shown in the diagrams correspond to the following models of the signal: f D = 500 
Hz, r = 40 yr, si = 4, s 2 = 3, s 3 = (solid), f a = 500 Hz, r = 1000 yr, si = 2, s 2 = 1, s 3 = (dashed), 
fo = 100 Hz, t == 40 yr, si = 4, s 2 = 2, s 3 = (dotted), / = 100 Hz, r = 1000 yr, si = 2, s 2 = 1, 
S3 = (dotted/dashed). The dimensionless amplitude of the waveform h Q = 1.1 x 10~ 23 for the models 
with f Q = 500 Hz and h a = 4.2 x 10^ 25 for the models with f a = 100 Hz. On the two top panels the lines 
which correspond to the models with the same frequency /„ coincide. 
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Figure 3: Cumulative distribution functions of the simulated rms errors of the signal parameters for 
various models of the signal in the case of all-sky searches. The observation time T = 7 days. The 
initial Hanford detector is assumed in the simulations. The neutron star parameters are the same as in 
Figure 1. The lines shown in the diagrams correspond to the following models of the signal: f a = 500 
Hz, t = 40 yr, si = 2, s 2 = 1, s 3 = (solid), f Q = 500 Hz, r = 1000 yr, si = 1, s 2 = 1, s 3 = (dashed), 
fo = 100 Hz, t = 40 yr, si = 2, s 2 = 1, s 3 = (dotted), / = 100 Hz, r = 1000 yr, si = 1, s 2 = 1, 
s 3 = (dotted/dashed). The dimensionless amplitude of the waveform h Q = 1.1 x 10~ 23 for the models 
with / = 500 Hz and h Q = 4.2 x 10~ 25 for the models with f Q = 100 Hz. On the two top panels the lines 
which correspond to the models with the same frequency So coincide. 
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Figure 4: Cumulative distribution functions of the simulated rms errors of the signal parameters for 
the individual detectors in the case of directed searches. The observation time T Q = 120 days. The 
neutron star parameters are the same as in Figure 1. The model of the signal's phase is described by 
s\ = 4, s 2 = 3, and s 3 — 0. The lines given on the diagrams correspond to various detectors: advanced 
Hanford (dotted/double dashed), initial Hanford (solid), VIRGO (dotted), wideband GEO600 (dashed), 
narrowband GEO600 (double dotted/dashed), and TAMA300 (dotted/dashed). 
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Figure 5: Cumulative distribution functions of the simulated rms errors of the signal parameters for 
directed searches and various models of the signal. The observation time T = 120 days. The initial 
Hanford detector is assumed in the simulations. The neutron star parameters are the same as in Figure 
1. The lines shown in the diagrams correspond to the following models of the parameter space: f D = 500 
Hz, r = 40 yr, si = 4, s 2 = 3, s 3 = (solid), f a = 500 Hz, r = 1000 yr, Sl = 2, s 2 = 1, s 3 = (dashed), 
f = 100 Hz, t == 40 yr, Sl = 4, s 2 = 2, s 3 = (dotted), /„ = 100 Hz, r = 1000 yr, si = 2, s 2 = 1, 
S3 = (dotted/dashed). The dimensionless amplitude of the waveform h a = 1.1 x 10 -23 for the models 
with f a = 500 Hz and h Q = 4.2 x 10~ 25 for the models with f a = 100 Hz. On the two top panels the lines 
which correspond to the models with the same frequency f a coincide. 
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Figure 6: Cumulative distribution functions of the simulated rms errors of the signal parameters for 
directed searches and various models of the signal. The observation time T = 7 days. The initial 
Hanford detector is assumed in the simulations. The neutron star parameters are the same as in Figure 
1. The lines shown in the diagrams correspond to the following models of the parameter space: f a = 500 
Hz, r = 40 yr, si = 2, s 2 = 1, s 3 = (solid), f Q = 500 Hz, r = 1000 yr, si = 1, s 2 = 1, s 3 = (dashed), 
fo = 100 Hz, t = 40 yr, si = 2, s 2 = 1, s 3 = (dotted), f a = 100 Hz, r = 1000 yr, si = 1, s 2 = 1, 
S3 = (dotted/dashed). The dimensionless amplitude of the waveform h a = 1.1 x 10 -23 for the models 
with / = 500 Hz and h Q = 4.2 x 10~ 25 for the models with f„ = 100 Hz. On the two top panels the lines 
which correspond to the models with the same frequency f a coincide. 
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Figure 7: Cumulative distribution functions of the simulated rms errors of the signal parameters for all- 
sky searches including the two proper motion parameters [i a and fig. The observation time T a = 120 days. 
The advanced Hanford detector was taken in the simulation. We assume that star's ellipticity e = 10~ 5 , its 
moment of inertia w.r.t. the rotation axis / = 10 g cm 2 , its distance from the Earth r Q = 40 pc, and the 
frequency f a = 500 Hz [these values give the dimensionless amplitude of the waveform h a = 2.6 x 10~ 22 ]. 
The spindown age r = 40 years. We have assumed that neutron star transverse velocity is 10 3 km/s. In 
the bottom right panel we give the distributions of the rms errors cr(/x Q ) (solid) and o-(ps) (dotted). 
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Figure 8: Comparison of cumulative distribution functions of the simulated rms errors of the signal 
parameters in the case of all-sky searches for the three simplified models of the gravitational- wave signal 
with the signal model given in Section 2. The observation time T = 7 days. The initial Hanford detector 
is assumed in the simulations. The neutron star parameters are the same as in Figure 1 so the solid lines 
are the same as solid lines in Figure 1. The other lines are as follows: the constant amplitude model 
(dotted), the linear model I (dotted/dashed), and the linear model II (dashed). 
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Figure 9: Comparison of cumulative distribution functions of the simulated rms errors of the signal 
parameters in the case of all-sky searches for the three simplified models of the gravitational- wave signal 
with signal model given in Section 2. The observation time T a — 120 days. The initial Hanford detector 
is assumed in the simulations. The neutron star parameters are the same as in Figure 1 so the solid lines 
are the same as solid lines in Figure 1. The other lines are as follows: the constant amplitude model 
(dotted), the linear model I (dotted/dashed), and the linear model II (dashed). 
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Figure 10: Dependence of the rms errors of the phase parameters on the observation time T a . We have 
used the constant amplitude linear phase signal defined in Eqs. (|57j) and (|58|). The signal-to- noise ratio 
d = 10. We have assumed the latitude A of the Hanford LIGO detector, we have also put <fi r = 1.456 
and 4> = 0.123. The lines (and the numbers) in the diagrams correspond to different numbers s of 
spindown parameters included in the phase (J58|) : s — 4 (solid), s = 3 (dotted), s = 2 (dashed), s = 1 
(dotted/dashed), and s = (double dotted/dashed). The horizontal dotted lines give the rms errors of 
the parameters as predicted by the polynomial phase model with s = 4 spindowns included (see Appendix 
C). The numbers corresponding to these lines are taken from Table 4 of Appendix C (for the signal-to- 
noise ratio d = 10). We see that for the observation times T a > 1000 days the horizontal lines coincide 
with the solid lines. 
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Figure 11: Dependence of the rms errors of the phase parameters on the initial phase 4> r of the Earth's 
diurnal motion. We have used the constant amplitude linear phase signal defined in Eqs. ( f37| ) and ([58]) 
with s = 4 spindown parameters included. The observation time T Q — 7 days and the signal-to-noise ratio 
d = 10. We have assumed the latitude A of the Hanford LIGO detector, we have also put cj) = 0.123. 
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Figure 12: Dependence of the rms errors of the phase parameters on the initial phase 4> a of the Earth's 
orbital motion. We have used the constant amplitude linear phase signal defined in Eqs. ( p7| ) and ([58]) 
with s — 4 spindown parameters included. The observation time T Q = 120 days and the signal-to-noise 
ratio d = 10. We have assumed the latitude A of the Hanford LIGO detector, we have also put <fr r = 1.456. 
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Figure 13: Dependence of the rms errors of the phase parameters on the latitude A of the detector's 
location. We have used the constant amplitude linear phase signal defined in Eqs. ( p7| ) and (|5|) with 
s = 4 spindown parameters included. The observation time T a = 120 days and the signal-to-noise ratio 
d = 10. We have put 4> r = 1.456 and 4> Q = 0.123. 
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Figure 14: The rms errors a((k) of the phase parameters £ = ($o, ojq, . . . , uoa) plotted as functions of the 
initial time xo for the signal-to- noise ratio d = 10. The curves correspond to models of the phase with 
increasing number of spindown parameters: dashed/double dotted curve — no spindown paremeter (this 
corresponds to a monochromatic signal with unknown constant initial phase), dashed/dotted curve — 1 
spindown parameter, dashed curve — 2 spindown parameters, dotted — 3 spindown parameters, solid — 4 
spindown parameters. All the curves shown here are symmetric w.r.t. Xq ~ 0.5. For some discrete values 
of x curves for (r + l)-parameter model intersect with r-parameter model curves. This means that in 
such cases adding r + 1 parameter does not change the rms error in rth parameter. 
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Figure 15: Cumulative distribution functions of the simulated rms errors of the signal parameters for 
all-sky searches. The observation time T Q — 120 days. The advanced Hanford detector is assumed in 
the simulations. The neutron star parameters are the same as in Figure 1. For these set of parameters 
the model of the signal's phase consistent with the 1/4 of a cycle criterion is obtained from Eq. ( |i"2| ) 
by setting s% = 4, s 2 — 3, and S3 = 0. This is the model number 5 in the plots given here. Models 1 
through 4 are defined as follows: model 4 — s% = 3, s 2 = 3, S3 = 0; model 3 — s\ = 2, S2 = 2, s 3 = 0; 
model 2 — si = 1, s 2 = 1, s 3 = 0; model 1 — s\ = 0, s 2 = 0, s 3 = 0. We observe a curve merging of the 
various models that corresponds exactly to the intersection points of the curves in Figure 14 for x a = 0.5. 
This shows that qualitative properties of the Fisher matrix for the polynomial model are preserved in the 
exact model. The five cumulative distributions of the amplitude parameters in the two top panels are 
indistinguishable. This confirms an effective decorrelation of the amplitude and the phase parameters. 
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